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1.  ABSTRACT 


Experiments  were  conducted  to  investigate  the  mechanism  responsible  for  vortex  breakdown 
over  delta  wing.  By  varying  the  free  stream  velocity,  unsteady  vortex  breakdown  occurs  over  the 
wing  due  to  the  alteration  of  the  balance  between  leading-edge  vorticity  production  and 
convection.  This  allows  the  study  of  vortex  breakdown  in  both  the  spatial  and  temporal  domain. 
It  was  found  that  parameters  such  as  swirl  angle  and  circulation  do  not  correlate  well  with 
breakdown  when  the  flow  is  unsteady.  However,  vortex  breakdown  seemed  to  follow  the 
positive  feedback  mechanism  proposed  by  Brown  and  Lxipez.  When  breakdown  occurs,  there  is 
a  redistribution  of  vorticity  in  the  vortex  subcore,  where  the  vorticity  vector  changes  from  an 
essentially  axial  direction  to  a  generally  tangential  direction.  Further,  the  helix  angle  for  the 
velocity  exceeds  that  of  the  vorticity  globally  in  the  vortex  subcore  only  in  the  instant  just  before 
breakdown.  Finally,  it  was  found  that  through  breakdown,  global  quantities  such  as  circulation 
and  flow  rate  varies  in  phase  as  the  upstream  vorticity  and  velocity  changes  with  time.  This 
implies  that  the  breakdown  process  is  essentially  an  inviscid  mechanism. 

2.  OBJECTIVE 

Since  the  observation  of  the  breakdown  of  the  leading-edge  vortices  over  a  delta  wing  by 
Peckham  and  Atkinson  in  1957,  the  phenomenon  of  vortex  breakdown  remains  a  challenging  and 
unsolved  problem  in  aerodynamics.  This  problem  is  especially  important  for  today  s  high-agility 
aircraft  operating  at  large  angle  of  attack,  since  vortex  breakdown  above  the  wing  results  in 
abrupt  change  of  the  aerodynamic  characteristics,  general  unsteadiness,  and  could  lead  to  loss  of 
control.  Despite  the  many  theoretical,  experimental,  and  computational  studies  in  the  past,  the 
fundamental  mechanism  responsible  for  the  breakdown  of  the  leading-edge  vortices  is  not  clearly 
understood.  Many  of  the  theoretical  models  fails  to  explain  the  physical  mechanism  that  brings 
about  vortex  breakdown.  Most  of  the  theories  can  predict  vortex  breakdown  to  occur  within  the 
observed  range  of  swirl  angles  and  the  sensitivity  of  breakdown  to  an  axial  pressure  gradient,  but 
none  of  the  theories  can  predict  the  flow  in  the  breakdown  region  and  the  location  of  breakdown 
with  an  accuracy  sufficient  to  be  checked  against  experimental  measurements.  However,  if  the 
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fundamental  mechanism  for  breakdown  can  be  clarified,  the  insight  could  lead  to  prediction  of 
breakdown  position  and  the  development  of  appropriate  control  techniques. 

The  unsteady  water  channel  at  UCLA  provides  an  unique  experimental  apparatus  to  study 
vortex  breakdown.  The  water  channel  can  provide  highly-controlled  large-amplitude  velocity 
variation,  thus  enables  us  to  study  both  the  temporal  and  the  spatial  evolution  of  the  velocity  and 
vorticity  field  within  the  breakdown  region.  Combined  with  LDA  measurements  and  flow 
visualization  studies,  we  have  investigated  various  vortex  breakdown  criteria  and  have  attempted 
to  clarify  the  underlying  physical  mechanism  behind  vortex  breakdown. 

3.  VORTICAL  STRUCTURES  OVER  DELTA  WINGS 

Vortical  structures  over  delta  wings  result  from  boundary  layer  separation  along  the  leading- 
edges.  The  separations  generate  shear  layers  that  lead  to  the  formation  of  two  primary  vortices 
above  the  wing.  These  vortices  are  important  because  they  produce  up  to  60%  of  the  total  lift  at 

high  angle  of  attack  (Wentz  and  Kohlman  1971). 

The  flowfield  about  the  leading-edge  vortex  is  nearly  axisymmetric  around  the  vortical  axis. 
Contrary  to  confined  vortices  generated  with  guidevane  in  vortex  tubes,  vortex  generated  by  a 
swept  delta  wing  is  formed  from  the  roll-up  of  the  vortex  sheet  which  is  continuously  shed  from 
the  leading-edge,  so  that  a  significant  variation  with  axial  distance  along  the  vortex  is  an  inherent 
aspect  of  the  vortex  structure.  The  vortex  in  a  tube  is  highly  axisymmetrical,  radially  confined 
and  normally  subject  to  only  slow  axial  variations.  The  flowfield  over  the  delta  wing  can  be 
characterized  by  its  two  primary  velocity  components,  the  axial  component  V^(r,x),  and  the  swirl 
or  azimuthal  component  Vefr.x).  The  radial  component  K  is  always  small  except  when 
breakdown  occurs. 

Figure  1  shows  the  typical  flow  profile  before  and  after  vortex  breakdown.  When  there  is  no 
vortex  breakdown,  the  axial  velocity  has  a  jet-like  profile,  with  the  centerline  velocity  that  is 
almost  twice  the  freestream  value.  However,  when  breakdown  occurs,  axial  velocity  changes  to  a 
wake-like  profile,  with  a  deficit  in  the  centerline  velocity.  Contrary  to  axial  velocity,  which 
shows  significant  differences  before  and  after  breakdown,  the  swirl  velocity  does  not  deviate 


2 


much  from  the  norm.  In  fact,  it  is  difficult  to  determine  when  breakdown  has  occurred  from  the 
swirl  velocity  alone. 

The  flowfield  inside  the  spiraling  vortices  can  be  divided  into  two  main  regions:  a  viscous 
subcore  and  a  outer  inviscid  rotational  core.  In  the  outer  part  of  the  flow,  the  vortex  motion 

approximates  to  irrotational  behavior  for  which  the  swirl  velocity  varies  as  Ve  -  -  -  On  the  other 

hand,  near  the  vortex  centerline,  viscous  effects  are  important,  and  the  swirl  velocity  behaves  like 
solid  body  rotation  for  which  Vg  ~  r  (figure  2,  where  r,  is  the  radius  of  the  vortex  subcore).  For 

an  ideal  vortex  that  is  axisymmetric,  the  swirl  velocity  profile  should  be  anti-symmetric  about  the 
vortex  center-line.  However,  due  to  presence  of  the  delta  wing,  the  vortex  is  skewed.  As  we  will 
see  later  on,  the  viscous  core  is  not  axisymmetric.  However,  to  simplify  our  calculations,  the 
flowfield  is  assumed  to  be  axisymmetric.  In  subsequent  studies  3-D  measurements  will  be 

conducted  to  verify  the  validity  of  the  axisymmetric  assumption. 

The  radius  of  the  viscous  subcore,  r,-,  can  be  derived  from  the  swirl  velocity.  The  viscous 
subcore  radius  is  defined  as  the  distance  from  the  vortex  centerline  to  the  point  of  maximum 
swirl  velocity,  Vecmax)  (figure  2).  The  viscous  subcore  is  the  region  where  most  of  the  streamwise 

vorticity  is  confined.  Up  to  95%  of  the  streamwise  vorticity  is  contained  inside  the  viscous 
subcore. 

The  swirling  motion  of  the  leading-edge  vortices  can  also  be  characterized  by  the 
distribution  of  circulation,  F  =  2jtrVe .  The  circulation  has  a  parabolic  behavior  near  the  axis  and 

tends  towards  a  constant  value  F^  in  the  outer  part  of  the  flowfield.  However,  since  vorticity  is 
continuously  being  fed  into  the  vortex  by  the  delta  wing  leading-edge,  F ,  will  vary  both  in  the 
axial  and  radial  direction., 

A  common  parameter  used  to  describe  the  swirling  motion  is  the  swirl  angle  6,  where 


The  swirl  angle  depends  on  the  angle  of  attack,  the  sweepback  of  the  wing,  and  the  yaw  angle.  In 
general,  the  swirl  angle  increases  with  the  angle  of  attack,  decreases  with  sweepback  angle,  and 
increases  with  yaw  (Werle  1971).  The  swirl  angle  is  an  important  parameter,  as  many  studies 
with  vortex  tubes  have  found  that  breakdown  generally  occurs  when  the  swirl  angle  reaches  a 
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critical  value  (typically  greater  than  45  degrees)  at  some  point.  Figure  3  shows  typical  variations 
of  swirl  angle  and  circulation  with  radial  distance  from  the  vortex  centerline.  The  swirl  angle  is 
zero  at  the  center  of  the  vortex  core  and  reaches  a  maximimi  at  the  edge  of  the  subcore,  where  the 
swirl  velocity  is  maximum. 

During  vortex  breakdown,  since  there  is  an  abrupt  change  in  the  structure  of  the  vortex  core 
with  rapid  expansion  of  core  size,  there  is  also  a  rearrangement  of  the  spatial  distribution  of 
vorticity.  By  assuming  axisymmetric  flow,  the  axial,  azimuthal,  and  radial  components  of 
vorticity  in  cylindrical  coordinates-i2j|.,  and  D-czn  be  computed  from  the  axial  and  swirl 
velocities: 


=— + — ^ 
r  dr 


Dq  - 


dV,  dK 


n=- 

dx 


dx  dr 

dK 
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dF 

is  generally  very  small.  Further,  since  Vr  is  normally  very  small  as  compared  to  — — ,  the 

dr 

dV 

swirl  or  azimuthal  vorticity,  £2^  can  be  approximated  by  eliminating  — - .  Figure  4  shows  the 

vorticity  components  calculated  from  the  velocity  proflles  in  figure  1 .  Since  the  swirl  velocity 
does  not  change  much  during  transition  to  breakdown,  axial  vorticity  remains  high  in  the  vortex 
centerline.  The  majority  of  the  axial  vorticity  is  contained  within  the  vortex  subcore.  On  the 
other  hand,  we  observe  that  in  the  transition  to  breakdown,  the  swirl  vorticity  component  changes 

dV 

sign  with  breakdown.  The  highest  values  of  £2  occur  where  — —  gradients  are  the  largest. 

^  dr 

dV 

When  there  is  no  vortex  breakdown,  the  axial  velocity  profile  is  jet-like  and  — —  is  always 

dr 

negative,  hence  positive  £2  .  However,  when  vortex  breakdown  occurs,  axial  velocity  profile  is 


wake-like  with  a  deficit  in  the  center  part,  and 


dV 

X 

dr 


is  positive  inside  the  vortex  core,  thus 
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negative  Q^.  The  significance  of  the  production  of  negative  swirl  vorticity  will  be  discussed 
further  in  sections  below. 


4.  VORTEX  BREAKDOWN  MECHANISM 


4.1  Vorticity  Transport 


The  vortex  breakdown  mechanism  may  be  simplified  and  explained  by  the  concept  of 
vorticity  transport  developed  by  Reynolds  and  Carr  (1985).  They  proposed  that  in  order  for  the 
leading-edge  vortices  to  remain  steady,  it  is  necessary  that  the  axial  convection  of  vorticity  along 
the  vortex  core  balances  the  influx  of  vorticity  generated  by  the  separation  of  the  leading-edge 
boundary  layer.  Vortex  breakdown  can  be  partially  explained  as  a  disturbance  of  the  vorticity 
balance  due  to  a  reduction  of  the  axial  convection  of  vorticity. 

A  characteristic  measure  of  the  balance  between  vorticity  influx  and  convection  is  the  swirl 
angle.  If  the  angle  of  attack  is  increased  the  swirl  velocity  increases  as  well.  Simultaneously, 
axial  convection  of  vorticity  decreases  because  of  the  imposed  strong  adverse  pressure  gradient  at 
the  trailing  edge  resulting  from  the  Kutta  condition.  Consequently,  the  axial  location  of 
breakdown  moves  closer  to  the  apex  with  increasing  angle  of  attack.  The  same  trend  was 
observed  with  smaller  sweepback  angles  and  increased  yaw.  Studies  with  vortex  tubes  have 
found  that  breakdown  occurs  when  the  swirl  component  of  the  flow  was  sufficiently  large, 
typically  when  the  swirl  angle  is  above  45  degrees  (i.e.,  when  the  swirl  velocity  exceeds  the  axial 
velocity).  In  terms  of  vorticity  balance,  a  large  swirl  angle  indicates  that  the  vorticity  generated 
at  by  the  leading-edge  separation  is  not  been  convected  effectively  due  to  the  decreased  axial 
velocity. 

We  can  gain  insight  into  the  mechanism  for  initiating  vortex  breakdown  by  examining 
vorticity  transport  for  the  simplest  case.  The  vorticity  transport  equation  is: 

—  =  (:!2.v)f  +  vV=/7  (5) 


ignoring  viscosity,  the  axisymmetric  vorticity  transport  equation  in  the  azimuthal  direction  is: 
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dx 
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For  the  simplest  case,  assume  that  Vr  and  Qr  are  small,  Vff=V^x,r),  14=V/rX  we  get 

(7) 


dt  "  a;c 


Equation  (7)  describes  the  temporal  development  for  the  redistribution  of  the  axial  vorticity 
into  the  swirl  vorticity  component  if  the  swirl  velocity  component  decreases  in  the  axial 
direction  x,  Qq  is  negative  and  will  decrease  with  time.  If  0,  the  swirl  vorticity  distribution 
will  induce  an  axial  velocity  component  opposite  to  the  main  flow  direction  due  to  the  Biot- 
Savart  law.  This  causes  a  strong  deceleration  of  the  axial  velocity  component. 


4.2  Positive  feedback  Mechanism 

The  positive  feedback  mechanism  for  breakdown  in  pipes  has  been  proposed.  A  positive 
axial  pressure  gradient  at  large  radial  distances  from  the  axis  of  the  slender  vortex  leads  to  a 
deceleration  of  the  axial  velocity  component.  Conservation  of  mass  requires  radial  outflow  and 
from  the  conservation  of  angular  momentum  and  the  vorticity  transport  equation,  a  redistribution 
of  the  axial  vorticity  component  into  the  azimuthal  vorticity  component  follows.  The  azimuthal 
vorticity  component  induces  an  additional  axial  velocity  component  against  the  main  flow 
direction,  resulting  in  an  even  stronger  deceleration  of  the  axial  velocity,  and  therefore 
enhancement  of  the  radial  velocity  component.  This  feedback  amplifies  the  redistribution  of  the 
axial  vorticity  into  the  azimuthal  vorticity.  This  process  may  lead  to  the  formation  of  a 
stagnation  point  on  the  axis  of  the  slender  vortex,  followed  by  a  region  of  recirculating  flow. 

Brown  and  Lopez  (1990)  were  the  first  to  note  the  positive  feedback  loop  for  the  initiation  of 
bubble-type  vortex  breakdown  in  a  pipe.  The  above  argument  can  also  be  made  for  the 
development  of  the  vortex  ^  lines:  If  the  flow  on  the  axis  of  a  slender  vortex  is  decelerated, 
conservation  of  mass  requires  radial  outflow.  The  vortex  lines,  which  are  straight  and  parallel 
upstream  of  the  point  of  deceleration,  are  stretched  and  tilted  in  the  region  of  radial  outflow.  The 
tilting  of  the  axial  vorticity  is  a  consequence  of  the  conservation  of  circulation  along  a  stream 
surface  in  steady,  inviscid  flow.  The  vorticity  vector  Q  acquires  a  component  in  the  azimuthal 
direction,  and  the  axial  flow  is  decelerated  due  to  induction.  Stretching  and  tilting  of  the  vortex 
lines  are  further  enhanced.  This  nonlinear  interaction  between  stretching  and  tilting  and 
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deceleration  of  the  axial  flow  may  be  promoted  by  a  positive  axial  pressure  gradient  at  large 


radial  distances.  The  interaction  may  cause  an  upstream  migration  of  the  stagnation  point 
leading  to  either  a  stable  vortex  or  a  bubble-  or  spiral-type  breakdown.  The  mechanism  is 
assumed  to  be  purely  inviscid  and  relies  upon  a  radial  outflow  or,  equivalently,  an  adverse 
pressure  gradient,  to  drive  the  breakdown  process. 


In  the  absence  of  viscous  or  turbulent  diffusion.  Brown  and  Lxipez  postulated  that  a 
necessary  condition  for  vortex  breakdown  to  occur  is  that  the  helix  angle  for  the  velocity  exceeds 

Ve 

that  of  the  vorticity  on  some  stream  surfaces.  That  is,  on  some  stream  surfaces,  y  ^  q  • 


4.3  Inviscid  or  Viscous  Mechanism 

In  general,  most  theoretical  studies  of  vortex  breakdown  may  be  classified  into  the  following 
categories: 

1.  Quasi-cylindrical  approach  and  analogy  to  two-dimensional  boundary  layer  separation 

(Hall  1965). 

2.  Hydrodynamic  instability  (Ludwieg,  1962). 

3.  Wave-motion  theories: 

a.  Finite  transition  between  critical  states,  analogous  to  a  hydraulic  jump  (Squire  1960, 
Benjamin  1962). 

b.  Breakdown  is  the  re.sult  of  the  trapping  of  long,  weakly  nonlinear  wave  propagating 
in  nearly  critical  swirling  flows  (Randall  and  Leibovich  1973). 

c.  Two-stage  transition  (Escudier  and  Keller  1 983). 

Flow  visualization  has  shown  that  the  structure  and  location  of  the  leading-edge  vortices  is 
independent  of  the  Reynolds  number  (Lee  and  Ho,  1990).  This  seems  to  imply  that  the 
breakdown  is  governed  by  an  inviscid  mechanism.  However,  the  detailed  flow  structure,  such  as 
the  location  and  strength  of  the  secondary  flow  separation,  depends  on  whether  the  leading  edge 
is  sharp  or  round.  Furthermore,  the  viscous  core  of  the  leading  edge  vortices  decreases  in  size 
with  increasing  Reynolds  number  due  to  the  thinning  of  the  boundary  layer  at  the  leading  edge. 
In  addition,  while  a  few  inviscid  theories  may  invoke  viscosity  to  initiate  breakdown,  all  inviscid 
theories  contend  that  the  rapid  transformation  of  the  vortex  cannot  be  attributed  to  viscosity.  It  is 
rather  a  result  of  inertial  effects  such  as  instability  or  wave  propagation. 
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On  the  other  hand,  viscous  theories  point  to  the  slenderness  of  the  vortex  subcore  and  argue 
that  viscous  forces  play  the  dominant  role  for  its  generation.  Since  the  subcore  is  slender, 
viscous  diffusion  could  occur  rapidly.  Using  the  analogy  to  two-dimensional  boundary  layer 
separation,  since  the  gradients  of  vorticity  becomes  increasingly  large  towards  the  vortex  axis,  it 
is  expected  that  viscous  diffusion  could  become  appreciable.  By  feeding  the  outer  mviscid 
solution  into  the  full  Navier-Stokes  equation  and  comparing  the  inertial  and  viscous  terms.  Hall 
(1965)  argues  that  inside  the  subcore  viscous  forces  are  as  important  as  inertial  forces.  The 
governing  equations  demonstrates  that  the  interplay  between  radial  inertial,  pressure,  and  viscous 
forces  are  important  inside  the  subcore  (Berger  and  Erlebacher  1995).  Consequently,  some 
argues  that  breakdown  is  necessarily  a  viscous  phenomenon. 

Experimental  results  that  can  shed  more  light  on  the  validity  of  these  theories  are  not 
available  as  the  flow  regime  is  three-dimensional  and  involves  large  flow  gradients  within  a 
small  region.  Hence  any  measurement  technique  has  to  be  not  only  non-intrusive  but  also 
sufficiently  accurate  to  resolve  the  important  flow  features.  Furthermore,  since  vortex 
breakdown  is  an  axially  evolving  flow  field,  measurements  would  need  to  be  taken  at  several 

axial  locations  with  accurate  spatial  resolution. 

The  above  problems  were  overcomed  in  the  present  work  on  a  delta  wing  in  unsteady  free 
stream.  In  this  novel  approach,  in  addition  to  the  spatial  domain,  the  mechanism  of  vorticity 
transport  over  a  delta  wing  was  analyzed  in  the  temporal  domain  by  varying  the  upstream 
production  of  leading-edge  vorticity.  This  enables  the  study  of  vortex  breakdown  both  in  the 
spatial  and  the  temporal  domain.  In  addition,  this  allows  the  breakdown  process  to  be  studied  in 
its  most  fundamental  definition-the  failure  of  the  fluid  to  transport  vorticity  in  an  organized  and 
efficient  fashion.  As  the  upstream  leading-edge  vorticity  underwent  large  amplitude  oscillations, 
the  vorticity  transport  within  the  core  changed  as  well,  even  to  the  point  of  breakdown.  The 
problem  thus  provided  the  means  for  ascertaining  whether  the  dominant  mechanism  of 
breakdown  was  viscous  or  inviscid.  If  viscous  diffusion  were  to  be  the  dominant  mechanism  of 
breakdown,  then  the  vorticity  transport  through  the  vortex  core  would  be  sharply  different  during 
and  after  breakdown.  On  the  other  hand,  if  the  breakdown  mechanism  was  inviscid,  there  will  be 
no  sharp  change  during  and  after  breakdown  as  the  vortex  core  would  still  continue  to  transport 
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vorticity  through  fluid  convection.  In  essence,  the  inviscid  mechanism  would  require  a  unique 
relation  between  the  fluid  flow  rate  and  the  fluid  circulation  within  the  vortex  core. 


5.  VORTEX  BREAKDOWN  IN  UNSTEADY  FREE  STREAM 

5.1  Unsteady  Flow  Field 

All  experiments  were  conducted  in  the  unsteady  water  channel  at  UCLA.  The  freestream 
velocity  varies  as: 

-^=l  +  i?coso)r  =  l  +  /?cos-^^  (8) 

t/„  T 

where  C/„  is  the  time-averaged  freestream  velocity,  R  is  the  dimensionless  amplitude  {R<\ ),  and 
(0=2ti/T\s  the  radial  frequency  of  the  velocity  variation.  The  reduced  frequency  can  be  defined 
as  k=coc/2U^,  and  is  the  ratio  between  the  vortex  convection  time-scale  (cAJJ  and  the  external 
perturbation  time  scale  (27t/0)). 

The  global  feature  of  the  vortex  breakdown  phenomenon  was  examine  by  flow  visualization 
for  both  the  steady  and  unsteady  cases.  Aspect  ratio  1,  2  and  4  (sweep  back  angle  76  deg,  63.4 
deg,  and  45  deg  respectively)  delta  wings  were  studied  at  different  angle  of  attack.  Vortex 
breakdown  depends  both  on  the  angle  of  attack  and  the  sweep  back  angle.  For  a  given  sweep 
angle,  vortex  breakdown  takes  place  near  the  trailing  edge  and  moves  upstream  with  increasing 
angle  of  attack  (note  steady  free  stream  pictures  in  figure  5).  At  the  same  angle  of  attack,  larger 
sweep  angles  delays  the  occurrence  of  breakdown  at  the  trailing  edge.  For  a  given  angle  of  attack 
and  sweep  angle,  when  steady  breakdown  position  is  near  the  midcore  and  the  trailing  edge, 
unsteady  forcing  can  cause  the  breakdown  position  to  vary  by  large  amplitudes,  as  large  as  40- 
50%  of  the  chord  length  during  an  unsteady  cycle  (cf.  figure  5  (a)  and  (b)  for  a=15  and  20 
degrees).  However,  for  very  small  and  very  large  angle  of  attacks,  i.e.,  when  the  steady 
breakdown  position  is  very  far  downstream  or  near  the  apex,  vortex  breakdown  locations  varies 
very  little  during  unsteady  forcing  (figure  5  (d)  and  (e),  a=30  and  40  degrees). 

For  the  cases  where  there  is  large  variation  of  the  vortex  breakdown  position  during 
unsteady  forcing,  flow  visualization  shows  that  the  breakdown  positions  moves  downstream 
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during  flow  deceleration  (0<t/T<05),  but  suddenly  appears  at  an  upstream  location  during 

flow  acceleration  (figure  5  (b),  A=2,  C35=20  degrees,  k=OJ66t  /?=0.583,  t/T'-O.lS),  and  then  moves 
downstream.  This  implies  that  the  axial  velocity  is  undergoing  a  change  in  flow  property. 

To  verify  this  transition,  detailed  velocity  field  measurement  were  conducted  along  a  plane 
across  the  vortex  core  for  an  aspect  ratio  one  delta  wing.  Measurements  of  velocity  components 
parallel  and  normal  to  the  wing  surface  were  taken  in  the  plane  containing  the  vortex  core  axis. 
Since  the  measurements  were  taken  along  a  traverse  line  parallel  to  the  trailing  edge,  the  velocity 
components  measured  approximately  correspond  to  the  axial  and  swirl  velocities  across  the 
vortex  core.  Phase-averaged  measurements  were  taken  at  different  streamwise  locations  for 
05=15  degrees,  amplitude  /?=0.42,  and  reduced  frequency  *=1.27.  The  origin  of  the  coordinate 
system  is  located  at  the  trailing  edge  and  centered  about  the  vortex  core.  Figure  6  shows  the 
spatial  variation  of  axial  velocity  through  an  unsteady  cycle.  When  breakdown  occurs  over  the 
wing,  the  axial  velocity  distribution  shows  an  abrupt  transition  from  a  jet-like  to  a  wake-like 
profile  at  the  breakdown  location.  During  flow  deceleration  (0«/r<0.5),  breakdown  moves 
downstream  beyond  the  trailing  edge  of  the  delta  wing.  After  the  breakdown  leaves  the  wing, 
axial  velocity  profile  everywhere  on  the  wing  is  jet-like  (t/T=0.375),  with  velocity  close  to  two 
times  the  freestream  velocity.  At  t/T=0.875,  axial  velocity  is  almost  uniform  everywhere  along 
the  wing.  Axial  velocity  profiles  then  transform  to  wake-like  at  most  axial  locations  after  this 
instant.  Figure  6  verifies  that  by  varying  the  freestream  velocity,  we  can  manipulate  the  vortex 
breakdown  position  as  to  enable  the  study  vortex  breakdown  not  only  in  the  spatial  domain,  but 
also  the  temporal  evolution  of  breakdown. 

As  we  have  noted  in  figure  1  ,  the  swirl  velocity  profile  does  not  change  much  through  the 
breakdown  process.  This  is  evident  in  the  swirl  velocity  surveys  (figure  7)  during  flow 
deceleration  and  acceleration.  The  swirl  velocity  profile  remains  similar  spatially  even  during 
transition  between  no-breakdown  and  breakdown.  Note  that  the  second  band  of  high  velocity 
region  beyond  r  ~  16  mm  is  already  beyond  the  outer  edge  of  the  delta  wing  leading-edge,  thus 
velocity  field  in  that  region  approximates  the  normalized  freestream  value.  Figures  8  is  the  same 
velocity  profiles  of  figures  6  and  7  plotted  in  vector  form. 
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The  swirl  velocity  can  be  used  to  determine  the  radius  of  the  vortex  viscous  subcore  (cf. 
Figure  2).  Figures  9  shows  the  variation  of  the  vortex  viscous  subcore  radius  during  the  unsteady 
cycle.  As  expected,  when  breakdown  occurs,  the  subcore  expands  in  size,  and  vice  versa. 

Studies  with  vortex  tubes  have  found  that  breakdown  occurs  when  the  swirl  component  of 
the  flow  was  sufficiently  large,  typically  when  the  swirl  angle  is  above  45  degrees  (i.e.,  when 
Contrary  to  those  studies,  the  present  data  shows  that  the  swirl  angle  is  not  a  good 
indicator  of  breakdown  when  the  flow  is  unsteady  (figure  10).  Comparing  axial  velocity 
distributions  with  swirl  angle,  there  is  no  direct  correlation  between  maximum  swirl  angle  and 
vortex  breakdown.  In  fact,  the  highest  swirl  angle  occurs  at  time  intervals  where  there  axial 
velocity  is  jet-like  everywhere  in  the  plane  (t/T=0.375). 

Since  dramatic  changes  take  place  in  the  vortex  core  during  the  breakdown  process,  its  effect 
on  global  quantities  such  as  circulation  is  of  primary  interest.  Ng  (1989)  proposed  that  the  vortex 
may  have  a  limit  on  the  maximum  amount  of  vorticity  per  unit  area  at  a  given  station. 
Consequently,  a  critical  vorticity  concentration  occurs  above  which  the  aerodynamic  forces 
cannot  maintain  a  stable  vortex  over  the  airfoil.  If  this  concentration  is  exceeded,  the  vortex 
transitions  to  another  state,  such  as  post  breakdown,  to  redistribute  the  excess  vorticity.  In  other 
words,  the  local  flow  condition  around  breakdown  is  similar  to  that  of  a  two-dimensional  airfoil 
during  static  stall.  For  example,  an  increase  in  the  angle  of  attack  leads  to  a  higher  rate  of 
generation  of  the  axial  vorticity  component  without  an  accompanying  increase  in  the  axial 
velocity.  The  subsequent  increase  in  the  vortex  strength,  therefore,  leads  to  vortex  breakdown. 
However,  figure  1 1  shows  that  the  circulation  is  not  influenced  by  the  breakdown,  as  the  profiles 
seems  to  be  similar  and  varies  with  freestream  velocity  variation. 

Axial  vorticity  calculations  (Figure  12)  shows  that  most  of  the  vorticity  is  contained  within 
the  vortex  viscous  subcore,  and  varies  very  little  with  breakdown.  This  is  expected  since  the 
swirl  velocity  does  not  give  indication  of  breakdown.  Therefore,  simply  determining  the 
maximum  value  of  the  axial  vorticity  does  not  make  a  definitive  statement  of  the  status  of  the 
vortex. 

The  Brown-Lopez  hypothesis  argues  that  the  breakdown  mechanism  relies  on  the  production 
of  negative  azimuthal  vorticity.  This  results  from  a  tilting  and  stretching  of  the  predominantly 
axial  vorticity  vector  £2x-  During  breakdown,  the  mean  rotational  vector  Q  changes  from  an 
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essentially  axial  direction  to  a  generally  tangential  direction  during  breakdown.  Figure  13  shows 
the  production  of  negative  azimuthal  vorticity  in  the  regions  where  breakdown  occurs  (cf.  figure 
6). 

To  better  examine  how  the  vorticity  is  being  redistributed  during  breakdown,  the  amount  of 
vorticity  has  been  computed  by  integrating  the  axial  and  azimuthal  vorticity  ( H2A  /  cU^ )  in  the 
inner  region  (-16  <  r  <  16  mm).  Since  most  of  the  vorticity  is  contained  in  the  viscous  subcore, 
and  the  subcore  radius  is  always  less  than  16  mm,  the  integrated  vorticity  values  will  give  us  an 
idea  of  the  total  amount  of  vorticity  flux.  Figure  14  shows  the  temporal  variation  of  the 
integrated  vorticity  at  different  streamwise  locations.  The  dash  lines  denote  negative  vorticity.  It 
can  be  clearly  seen  that  when  breakdown  occurs,  the  vorticity  has  been  distributed  from  the  axial 
component  to  the  production  of  azimuthal  vorticity.  Furthermore,  at  the  upstream  position  where 
breakdown  does  not  occur  (X  =  -100  mm),  the  axial  vorticity  varies  as  the  freestream  velocity. 
However,  at  the  downstream  position  (X  =  0  mm),  there  is  a  phase  shift  of  the  integrated  vorticity 
values  with  freestream  velocity.  This  is  due  to  the  delay  caused  by  the  vortex  convection  from 
the  leading-edge.  This  will  be  discussed  in  more  detain  in  the  section  below.  In  addition,  it  is 
difficult  to  correlate  between  breakdown  and  the  maximum  amount  of  vorticity,  as  the  present 
data  does  not  provide  evidence  to  support  the  argument  that  breakdown  occurs  as  a  result  of  a 
total  amount  of  vorticity  being  reached  at  a  given  streamwise  position. 

In  addition  to  the  rotation  of  the  vorticity  vector  and  the  production  of  negative  azimuthal 
vorticity.  Brown  and  Lopez  postulated  that  a  necessary  condition  for  vortex  breakdown  to  occur 
is  when  the  helix  angle  for  the  velocity  exceeds  that  of  the  vorticity  on  some  stream  surfaces,  i.e.. 
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Figure  15  shows  the  Brown  and  Lopez  criterion  at  different  phases  of  the  unsteady 


v;  Qf. 

cycle.  The  criterion  has  been  dichotomized  such  that  “>7^  is  shown  in  white  and 


—  <  — ^  is  shown  in  black. 


Comparing  with  the  radius  of  the  vortex  subcore  (figure  9),  we 


find  that  the  Brown-Lopez  criterion  is  globally  satisfied  inside  the  vortex  subcore  only  just  before 
breakdown  (t/T=0.875).  The  criterion  holds  true  at  all  radial  locations  inside  the  subcore  only  at 
j/r=0.875  and  not  at  any  other  instant.  In  the  outer  flow  region,  the  gradients  of  axial  and  swirl 
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velocities  are  small,  thus  susceptible  to  large  errors  when  computing  and  comparing  vorticity 
ratios.  Consequently  for  the  present  case  we  can  only  verify  the  Brown-Lopez  condition  inside 
the  subcore,  where  the  velocity  gradients  are  large.  Nevertheless,  our  present  data  does  verify  the 
validity  of  the  Brown-Lopez  condition  inside  the  vortex  subcore. 

5.2  Unsteady  Vortex  Breakdown  Mechanism 

The  organization  of  the  unsteady  flowfield  can  also  be  revealed  by  examining  the 
streamlines,  which  are  defined  as  iso-value  lines  of  the  stream  function  *F(X,r),  where, 

nx,r)=^^^dr\  (9) 

Further,  since  the  flowrate  for  axisymmetric  flow  is  defined  as  Qm(X.r),  where 

|2.(X,r)  =  2^|;-^T,dn  (10) 

the  streamlines  also  represent  the  constant-mass  flow  lines.  Figure  16  shows  the  unsteady  flow 
streamlines.  The  constant-mass  flow  lines  have  been  normalized  by  the  maximum  mass-flow  rate 
at  the  outer  flow.  Each  successive  streamline  outwards  from  the  vortex  center  (r=0)  has  twice 
the  mass-flow  rate  as  the  previous  one.  Throughout  the  unsteady  cycle  the  streamlines  stay 
straight  and  parallel  in  the  axial  direction,  verifying  that  the  vortex  flow  is  quasi-cylindrical. 
When  there  is  breakdown  over  the  delta  wing  (the  duration  from  t/T=0.875  to  1 .250),  there  is 
slight  divergence  of  the  streamlines  along  the  axial  locations  where  breakdown  is  occurring. 
Conservation  of  mass  requires  this  radial  outflow  due  to  axial  velocity  deceleration.  To  see  how 
the  flow  rate  is  evolving  in  the  time  domain,  figures  1 7  examine  how  the  flow  rate  changes  at 
each  axial  location  during  the  unsteady  cycle.  Again,  at  the  axial  locations  where  breakdown 
occurs,  the  constant-mass  flow  rate  contract  and  expands  with  upstream  flow  condition.  Figure 
1 8  shows  the  temporal  variation  of  vortex  viscous  core  size  (r,.)  at  different  axial  locations.  The 
viscous  subcore  size  expands  and  contract  in  temporally  in  response  to  breakdown  and  upstream 
velocity  variations. 

The  temporal  evolution  of  circulation  (figure  19)  also  shows  that  the  lines  of  constant 
circulation  at  any  give  location  contract  and  expand  rhythmically.  To  get  a  better  picture  of  how 
the  circulation  at  a  fixed  radial  position  changes  with  time,  the  figure  1 9  has  been  replotted  as 
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surface  plots  at  selected  axial  locations  (figure  20).  The  circulation  at  any  radial  distance  from 
the  vortex  core  exhibits  cyclical  variation  with  time.  Beyond  r  ~  20  mm,  the  measurements  are 
already  beyond  the  outer  leading-edge  of  the  delta  wing,  thus  the  circulation  varies  in  phase  with 
the  freestream  velocity.  However,  as  we  move  inward  towards  the  vortex  core,  the  variation  of 
circulation  has  a  phase-lag  with  the  free  stream.  This  delay  is  clearly  a  result  of  the  time  it  takes 
for  the  upstream  vorticity  to  convect  inside  the  core.  Vorticity  balance  concept  dictates  two  time 
scales  when  the  freestream  velocity  is  unsteady.  First,  the  change  in  the  vorticity  generation 
along  the  leading-edge  due  to  the  unsteady  freestream  is  transmitted  into  the  core  in  one  local 
vortex  turn-over  time  (the  characteristic  time  for  one  local  revolution  of  the  vortex).  For  our 
case,  we  can  say  that  this  turn-over  time  is  closely  related  to  our  flow  forcing  time  scale  (27t/co). 
This  time  is  minimal  at  the  apex,  and  increases  towards  the  trailing  edge.  Therefore,  the  flow 
around  the  apex  is  more  sensitive  to  disturbances  than  at  any  other  location  along  the  leading 
edge.  This  explains  that  why  flow  visualization  shows  that  at  high  angle  of  attack,  when  the 
breakdown  is  near  the  apex,  often  asymmetric  vortices  are  observed  over  the  wing  (figure  5  (d)). 
The  second  time  scale  is  determined  by  the  time  during  which  an  upstream  disturbance  is 
convected  throughout  the  vortex.  This  characteristic  time  is  c/U^  where  c  is  the  mean  chord  of 
the  delta  wing  and  U„  is  the  average  freestream  velocity.  This  relationship  can  be  represented  by 
the  reduced  frequency.  Since,  in  general,  the  vortex  convection  time  scale  is  longer  than  the 
local  turn-over  time  (for  our  case  Jt=1.27),  streamwise  convection  becomes  the  limiting  factor  for 
the  vortex  to  respond  to  any  disturbance.  Because  of  the  streamwise  convection  time,  a  phase  lag 
of  response  occurs  relative  to  the  freestream  velocity  variation.  This  is  what  we  observe  in  figure 
20. 

It  is  interesting  to  note  that  the  viscous  core  size  expands  and  contracts  by  rather  large 
amplitude  with  time  (figure'  18).  However,  despite  the  increase  in  the  vortex  core  size,  the  flow 
rate  and  circulation  do  not  show  corresponding  large  variation.  K  viscous  diffusion  was 
responsible  for  transporting  the  vorticity  outwards,  we  should  see  circulation  and  flow  rate 
change  as  well.  This  seem  to  imply  that  most  of  vorticity  is  being  convected  away.  Further, 
although  the  flow  rate  does  not  show  very  big  change,  figures  17  and  19  show  that  the  flow  rate 
and  circulation  vary  in  phase  as  the  upstream  vorticity  and  velocity  changes  with  time.  This 
demonstrates  that  vorticity  is  convecting  with  the  fluid  particles  both  before  and  after 
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breakdown.  If  viscous  diffusion  were  to  have  played  a  dominant  role  during  breakdown,  then  the 
circulation  profile  should  show  a  difference  from  the  profile  for  flow  rate.  The  absence  of  any 
such  deviation  leads  to  the  conclusion  that  vortex  breakdown  is  an  inviscid  mechanism. 

6.  FUTURE  DIRECTIONS 

Previous  studies  of  vortex  breakdown  have  tried  to  formulate  and  quantify  flow  properties  as 
to  correlate  with  breakdown.  It  was  generally  thought  that  breakdown  occurs  when  certain 
parameters  reaches  some  critical  value,  such  as  the  swirl  angle  or  the  amount  of  vorticity.  As  we 
have  shown  for  unsteady  flows,  vortex  breakdown  is  an  extremely  complex  phenomenon  and 
cannot  be  characterized  by  any  single  parameter.  Vortex  breakdown  mechanism  involves  the 
local  properties  of  the  fluid  velocity  distributions,  as  well  as  global  properties  such  as  circulation 
and  flow  rate. 

Although  much  can  be  understood  from  a  single  cut  through  the  vortex  core,  more  work  is 
needed  to  evaluate  the  three-dimensional  flowfield  to  examine  the  global  influence  of  vortex 
breakdown.  In  the  present  studies  we  have  assumed  that  the  radial  components  of  velocity  and 
vorticity  is  small.  However,  that  assumption  might  not  be  valid  when  breakdown  occurs.  To 
really  clarify  the  role  of  viscous  diffusion  inside  the  vortex  subcore,  full  three-dimension  velocity 
and  vorticity  fields  will  have  to  be  considered.  The  Brown-Lopez  criterion  relies  on  the 
assumption  that  the  initial  reduction  in  the  azimuthal  component  of  vorticity  to  be  due  to  viscous 
diffusion  and  the  resultant  stretching  and  tilting  of  axial  vorticity.  As  the  flow  begins  to  diverge, 
the  further  production  of  negative  azimuthal  vorticity  is  dominated  by  inviscid  mechanisms. 
Detailed  three  dimensional  data  inside  the  vortex  core  will  enable  us  to  solidify  the  breakdown 
mechanism. 
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Fig  1.  Normalized  velocity  profiles  before  (t/T=0.375),  transition 
to  (t/T=0.875),  and  after  (t/T=l  .000)  vortex  breakdown. 
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Fig.  2.  Swirl  velocity  distribution 
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Fig.  3.  Variation  of  flow  property  distribution  with  radial  distance  from  the  vortex  centerline 
before  (t/T=0.375),  transition  to  (t/T=0.875),  and  after  (t/T=l  .000)  vortex  breakdown. 
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before  (t/T=0.375),  transition  to  (t/T=0.875),  and  after  (tfT=1.000)  vortex  breakdown. 
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Fig.  5  (a).  Flow  visualization  for  A  =  2  delta  wing,  a  =  15  deg,  k  =  0.837,  R  =  0.593. 
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Fig.  5  (b).  Flow  visualization  for  A  =  2  delta  wing,  a  =  20  deg,  k  =  0.766,  R  =  0.583. 
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Fig.  5  (c).  Flow  visualization  for  A  =  2  delta 
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Fig.  5  (d).  Flow  visualization  for  A  =  2  delta  wing,  a  =  30  deg,  k  =  0.822,  R  =  0.589. 
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Fig.  5  (e).  Flow  visualization  for  A  =  2  delta  wing,  a  =  40  deg,  k  =  0.8 1 1 ,  R  =  0.58 1 . 
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Fig.  6  (a).  Normalized  axial  velocity  survey  of  a  plane  across  the  vortex  core 
centerline  during  flow  deceleration.  a=15  deg,  A  =  1,  k  =  1.27,  R  =  0.42 
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Fig.  7  (b).  Normalized  swirl  velocity  survey  of  a  plane  across  the  vortex  core 
centerline  during  flow  acceleration.  (X=15  deg,  A  =  1,  k  =  1.27,  R  =  0.42 
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Fig.  8.  Velocity  vectors  in  a  plane  across  the  vortex  core  centerline  during  an 
unsteady  cycle:  (a)  normalized  axial  velocity,  (b)  normalized  swirl  velociy. 
(X?=15  deg.,  A  =  1,  k=1.27,  R  =  0.42. 
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Fig.  8,  Con’t.  Velocity  vectors  in  a  plane  across  the  vortex  core  centerline  during 
an  unsteady  cycle:  (a)  normalized  axial  velocity,  (b)  normalized  swirl  velociy. 
(X=15  deg.,  A  =  1,  k=1.27,  R  =  0.42. 
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Fig.  8,  Con’t.  Velocity  vectors  in  a  plane  across  the  vortex  core  centerline  during 
an  unsteady  cycle:  (a)  normalized  axial  velocity,  (b)  normalized  swirl  velociy. 
0=15  deg.,  A  =  1,  k=1.27,  R  =  0.42. 
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Fig.  9  (b).  Variation  of  vortex  viscous  core  position  during  acceleration. 
a=15  deg,  A  =  1,  k  =  1.27,  R  =  0.42 
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Fig.  10  (a).  Swirl  angle  survey  of  a  plane  across  the  vortex  core  centerline 
during  flow  deceleration.  a=15  deg,  A  =  1,  k  =  1.27,  R  =  0.42 
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Fig.  10  (b).  Swirl  angle  survey  of  a  plane  across  the  vortex  core  centerline 
during  flow  acceleration.  ot=15  deg,  A  =  l,k=1.27,R  =  0.42 
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Fig.  1 1  (b).  Normalized  circulation  survey  of  a  plane  across  the  vortex  core 
centerline  during  flow  acceleration.  a=15  deg,  A  =  l,k  =  1.27,  R  =  0.42 


Normalized  axial 
vorticity 


rVe  SVe'l  c 


leading-edge 


V  r 


0r 


U, 


-1 


-4a30-2ai00  102030 


trr  =  0.375 


E 

E 

X 


-40-3a2ai0  0  102030 


-40-30-20-10  0  102030 


t/T  =  0.125 


t/T  =  0.500 


-40-30-20-10  0  102030 
r  (mm) 


E 

E 

X 


-4a30-20-10  0  102030 
r  (mm) 


Fig.  12  (a).  Normalized  axial  vorticity  survey  of  a  plane  across  the  vortex  core 
centerline  during  flow  deceleration.  oc=15  deg,  A  =  1,  k  =  1.27,  R  =  0.42 
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Fig.  13  (a).  Normalized  azimuthal  vorticity  survey  of  a  plane  across  the  vortex 
core  centerline  during  flow  deceleration.  cx=15  deg,  A  =  1,  k  =  1.27,  R  =  0.42 
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Fig.  13  (b).  Normalized  azimuthal  vorticity  survey  of  a  plane  across  the  vortex 
core  centerline  during  flow  acceleration.  a=15  deg,  A  =  1,  k  =  1.27,  R  =  0.42 
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Fig.  14.  Temporal  variation  of  integrated  vorticity  values  at  different 
streamwise  locations.  a=  15  deg,  A  =  1,  k=1.27,  R  =  0.42. 
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Fig.  15  (a).  Brown  and  Lopez  vortex  breakdown  criterion  across  the  vortex 
centerline  during  deceleration.  (X=15  deg,  A  =  1,  k  =  1.27,  R  =  0.42 
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Fig.  15  (b).  Brown  and  Lopez  vortex  breakdown  criterion  across  the  vortex 
centerline  during  deceleration.  a=15  deg,  A  =  1,  k  =  1.27,  R  =  0.42 
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Fig.  17.  Temporal  variation  of  flowrate  across  the  vortex  core  centerline  at 
different  axial  locations.  a=  15  deg,  A  =  1,  k  =1.27,  R  =  0.42. 
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Fig.  17,  Con't.  Temporal  variation  of  flowrate  across  the  vortex  core  centerline 
at  different  axial  locations.  a=  15  deg,  A  =  1,  k  =1 .27,  R  =  0.42. 


ilf  Him 


0.00  0.20  0.40  0.60  0.80  1.00 

Time  (t/T) 

0.00  0.20  0.40  0.60  0.80  1.00 

Time  (t/T) 

X  =  -70  mm 

X  =  -100  mm 

0.00  0.20  0.40  0.60  0.80 

Time(t/T) 

X  =  -80  mm 


0.00  0.20  0.40  0.60  0.80  1.00 

Time(t/T) 


iO-f- 

0.00 


0.40  0.60 

Time(t/T) 


Fig.  19,  Con't.  Temporal  variation  of  normalized  circulation  across  the  vortex  core 
centerline  at  different  axial  locations.  cx=  15  deg,  A  =  1,  k  =1.27,  R  =  0.42. 


